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PREFACE 


Certain volumes of the TRACE documentation series were published by 
The Aerospace Corporation as Technical Operating Reports. Volume III: 
Trajectory Generation Equations and Methods was published as 
TOR-0066(9320)-2, Vol. IlI, and Volume V: Pifferential Correction “ro- 
cedure and Techniques was published as TOR-€66(9320)-2, Vol. V. 


Volume I: General Program Objectives, Description, and Summary was 
published as TR-0059(9320)-1, Vol. I, and Volume X: Lunar Gravity 
Analysis was published as TR-0059(9320)-1, Vol. X. Future volumes in 


this series will be published as Technical Reports. 
This report is published in two parts, A and P. 


The TRACE Program could not have cveen developed to its present status 
without the assistance of many xeople working in the fields of astrodynamics 
and software design. The authors acknowledge with gratitude and analysis 
and/or programming efforts of A. B. Bierman, R. J. Farrar, W. A. Feess, 
E. H. Fletcher, R. B. Freund, T. P. Gabbard, ©. G. Gibson, P. T. Gray, 
P. T. Guttman, J. A. Pearson, C. M. Price, W. F. Rearick, N. W. Rhodus, 
A. J. Rusick, L. J. Tedeschi, L. Wong, 1nd K. R. Young. In addition, 
consultations with W. T. Kyner, A. Troes:h, and H. H. Wertz have led to 


many significant improvements and added :.apabilities in the program. 


This report supersedes the appendices to TR-0059(9320)-1, Vol. VII, 
30 June 1971. WA at, BS 
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APPENDIX A 


DESCRIPTION OF THF GENERAL-PURPOSE 
INPUT ROUTINE GAIL1 


A.1 FUNCTION AND CAPABILITIES 


A.1.1 Purpose and Description 
* 
The purpose of GAIL1 is to read a set of Hollerith punched data and/or 


header cards into core with one FORTRAN cail statement. To achieve 

this, the data fields are converted to binary, while the header informa- 

tion is stored in Hollerith code; each piece of the datais converted and 
stored according to a specific conversion code. These data may be integers, 
real variables, arrays, alphanumeric data (short or long headers), ma- 


trices, or end-of-case indications. 


The data cards can be printed as received, or they can be placed ona 
scratch I/O device and then listed when an end of file (EOF) is detected on 
the input file. 


The subrovtine will always return control to the driver program except in 
the case of tape errors, system malfunctions, or completely erroneous 
punched data cards that disobey all FORTRAN I/O format rules. There- 
fore, the user must supply and test a status cell after an exit from the 


subroutine to determine which of the following occurred: 


e Normal end of case 
@ Some type of error 
e End of file reached on the input device 


“GAILI was originally developed by R. B. Gladson of The Aerospace 
Coiporation for the CDC 6000 computer series. It has been modified 
for IBM 360 and 370 computer series, 


A-le @ 


A.1.2 Method 


All data fields are converted according to FORTRAN I/O format 
specifications. Real variables are converted by the F specification, 
whether or not the exponent (E or D) is included with the variable. Inte- 
gers are converted by the I format, and alphanumeric data by the A format. 
This allows for alternate rnethcods of entering real variables and integers 
from the FORTRAN standards. This variety of metiiods (Sec. A. 1.3.7) 
enables the user to enter numbers in a manner more reasonable and 


mathematical than did previous input routines. 
A.1.3 Usage 
A.1.3.1 Data Card Format 


The data card format consists of three fields, each containing three sub- 
fields; a conversion code; a location; and a value (as in the following 
example): 


12 7 2728 33 5354 59 


S| es] ~ Als] 1b 


Data Field Data Field Data Field 


Subfield 1 2 3 
Conversion Code (C;) 1 27 53 
Location (L.) 2-6 28-32 54-58 
Value (V.) 7-26 33-52 59-78 


Columns 79 and 80 are not processed by the routine and can be used for card 
sequencing if the user desires. The conversion code (C,) is one of the 
alphabetic characters (Sec. A. 1.3.7); it specifies the type of conversion 

to ve used on the value field. The location subfield specifies the cell in 


which the converted value subfield is to be stored. 


‘ 


Aste 3:2 Header Card Format 


The header card forrnat consists of a conversion code in Column 1, a 
sequence number in Columns 2 through 6, anu any Hollerith information 


in Columns 7 through 78. Columns 79 and 80 are not used. 


t 2 ie 79 80 
ae ee a ee ee | 
Cc L Header seq 


A ts3 33 Decimal Points 


The use of a decimal point is valid in both integer and real variable 
entries; however, it is not mandatory to enter one. If the decimal point 
is missing from an entry, it is assumed to be a. the right of the last digit 


punched in the value field; for example: 


321. = 321 = 321E0 


A. 1.3.4 Signs 
If signs (+ or -) are used, they are placed in a column by themselves and 


immediately precede the value. Overpunches are not valid; e.g.: 


Valid 


Invalid 


a acai, DEPRES NS TAG A al 


eahiiiadinmaatetettine tabataciar gta ont Re re 


a 


i) 


Av1.355 Values (V. Subfield) 


It is not necessary for the entire subfield to be filled, since the first blank 
following the first nonblank denotes the end of the value. Superfluous low- 


order zeros should be omitted. 


If alphany meric -nformation is entered, it is stored in core as it appears 
in the value subdfiead. Thus, for this type of data entry, leading blanks 
are retained. If they appear, the blanks themselves are also stored in 


core, 


A. 1.3.6 Location (L. Subfield) 


The location is specified by either a variable, an array name, or the 


element subscripts in a single-dimension array or matrix. 


A.1.3.6.1 Variable, Array, or Matrix Name 


If the location contains at least one nonnumeric character that is not a 
comma as the third character, it is interpreted as a variable, array, or 
matrix name. The contents of the value subfield are stored in the cell for 
the variable or in the first cell of the array or matrix. This location then 
becomes the origin of all numeric and matr.x element locations that follow 
until another variable, array, or matrix name is encountered. Care must 


be taken to enter an array or matrix name before entering numeric or 


matrix element locations. 


A. 1.3.6.2 Single Numeric Entry 


If a location contains a single numeric value, the entry is treated as a 
subscript. With its previously established origin, this subscript is used 
to compute the location in which the entry in the value subfield is stored. 
For instance, if the established origin is NAME and the value in the loca- 
tion subfield of the input is 6, the entry in the value subfield is stored in 
cell NAME(6). (Single numeric entries must be left-justified, and leading 


zeros are optional.) In the following example, 


YT 2 7 
27 928 33 
53154 59 


LOCATION VALUE 


|_INAME 
Pale 


legitimate numeric entries wouid be <nown as: 


LOCATION 


- Not a valid subscript. 
- Store as 30, 


16) 
Qi - Not recognized as a relative 


address; it is flagged as an error, 


Sgt 


The symbol (1) indicates 1 space between digits, and @) indicates n spaces. 


A. 1.3.6.3 Matrix Element Entry 


A matrix element is stored by specifying two separate two-digit numeric 


entries, separated by a comma, in the location subfield. For example: 


-3.17 goes to element 
(3,4) of the matrix 


| +7.38 goes to element 
(15,1) of the matrix 


The elements are stored in the column major sort by using the formula 
@RIGIN (( - 1)1+ i) 


where I is the total number of rows in the matrix and (i, j) is the eleruent 


in question. 


A.1.3.6.4 Blank Location 


ee 


If the location is left blank, the entry in the value subfield is stored in the 
cell immediately following the one in which the last value was stored. 
Thus, an entire array may be entered by using a name in the location sub- 
field for the initial entry and leaving location subfields blank for ‘.e sub- 


sequent value subfield entries. The following is an exainple: 


LOCATION 


XYZ | stored in XYZ(1) 
— stored in XYZ(2) 


stored in XYZ (3) 


Ae 1 37 Conversion Code (Cc. £ubfield) 


The entry in the conversion code subfield specifies the method of converting 
the data in the value subfield. Care should be taken to apply the proper 
code with each item of data; if the proper code is not applied, a premature 
termination will occur ina FORTRAN format statement. The eight con- 


version codes currently available are described in the following subsections. 
A.1.3.7.1 Blank: Real Variable 
The number in the value subfield is converted ta and stored as a real 


variable, Entries can be made to the value subfield by any of four methods: 


Three of them are in the normal FORTRAN E, D, and F types of real 


variable input forrnats, e.g.: 


The fourth method allows for the input of real variables in scientific 
nomenclature (a X 10°), The number is written in the value subfield and 
is followed by an X, then by a 10, and then bv the power of ten to which 
the number is to be scaled (a sign must precede the exponent). Some 


examples follow: 


@ 14200 = 1. 4200x104 is input as 


or as 


@  .00762 = .762x107* is input as 


W) 


Aw 1.3. Tat 1: Integer 


The number in the value subfield is! converted to a fixed-point integer. A 
decimal point is allowable; if it is used, onlv the integer portion is con- 
verted. Before the value is converted by the FORTRAN I format, the 
decimal point and any digit(s) to the right are stripped off; e.g.: 


is stored as 71. 


A.1.3.7.3 D: Short Header (20 Characters) 


The contents of the value subfield are interpreted as 20 alphanumeric 
characters. They are stored in two consecutive cells in ascending order, 
beginning at the location specified by the location subfield. The charac- 
ters are stored in the same relative positions as in the value subfield; 


thus, leading, trailing, and embedded blanks are legitimate. For example: 


VALUE 


ABC(Q1)DEF(2)GHIJ(@) 234 


is stored in two consecutive words, with ABCQ)DEF@G in the first word, 
and HIJ@2 34 in the second. 


A. 1.3.7.4 H: Long Header (72 Characters) 


A card with an H in Column 1 is considered a long header card. If the 
location subfield is blank, the card is ignored (treated as a comments 
card). If the location subfield contains a left-justified integer 1, Columns 


7 through 78 of the card are stored directly in seven consecutive words in 


ascending core storage. 


Bal bi sR lita mS. Reet Wthaont 


The H conversion code is the only one for which a numeric location 
subfield is mandatory. If the location subfield is not blank or the integer 
{, the header is not stored, an error comment is printed, and the error 


termination procedure is invoked. 


A.1.3.7.5 L: Sparse Matrix Input Definitioa 


This input specifies only that matrix elements may be following, either 
on the same and/or on subsequent cards. It does not store any items of 
data; it only warns the subroutine that matrix elements might be stored 
and indicates where. The location subfield must contain the name of the 


matrix in which the items are to be located in core. 


The value subfield must contain two-digit integers separated by a comma. 

The entries are the size of the matrix [I and J ], and the comma 
max max 

must be the third character of the subfield. If it is not, an error comment 


is printed, and an error flag set. 


The Ls ie 


scripted locations of the matrix elements until they are redefined. Blank 


and J ase values are retained to compute the successive sub- 


location subfields may follow this definition if successive elements of the 
matrix are to be loaded. It must be remembered that the matrix is stored 


in column major. For example: 


VALUE 
’ 


indicates that a 3 xX 7 matrix G@OD is being defined and that following data 
values are to be stored in GOOD until another symbolic location name is 


encountered. Storage is to be made according to GOPD ((j - 1)I1+ i). 


A.1,3.7.6 M: Full Matrix Input Definition 


or: 


The use of this conversion code is identical to that of L, except that the 


entire matrix is preset to zero before anv of the elements are loadec. 
For example: 


ee se casio es 
Ys « 


33 so 
|ejrocation | vauve 
IMIG@@D_[03, 07 


a Necks set wit Saw 


presets the 3 x 7 matrix GO@D to zero; it is prepared to store the matrix 
elements foilowing it until another symbolic location name is encountered. 


4 A.1.3.7.7 Bs Octal 


The value subfieid is converted as a logical word. It is not necessary to 
include leading zeros, but the first octal digit must always occupy the left- 
QO most position in the subfield. If fewer than 20 digits are input, the value 


will be right- justified in the memory location. 


For example: 


| ¢| LOCATION VALUE 


is stored ac 00000000000000000123, in core. 


A.1.3.7.8 E: End of Case 


This defines the end of case, when control is returned to the object pro- 


gram from which it was called. The rest of the subfield and the remaining 
fields on the card are ignored. 


A.2 ERROR COMMENTS 


The input errors that the subroutine can recognize are listed in Table A-1t. 
In each case, an appropriate error comment is printed, the remaining 
cards in the data set are scanned until an E conversion code is found, and 


the run is terminated. The errors are not listed in any specific order. 


* 
. 


Table A-1. Errors Recognizable by Input Subroutine 


Error Message 


HEADER CARD NO. EXCEEDS MAXIMUM NO, 
ALLOWABLE, 


ERROR ATTEMPTING TO COMPUTE NEXT 
STORAGE LOCATION FROM BLANK LOCATION 
FIELD. BASIC SYMBOLIC ORIGIN HAS NOT 
BEEN SET. 


~— Field in Error 


ERROR ATTEMPTING TO COMPUTE NEXT 
STORAGE LOCATION FOR A MATRIX ELEMENT, 
BASIC SYMF DLIC ORIGIN HAS NOT BEEN SET. 


—— Field in Error 


MATRIX SIZE UNDEFINED. 


-— Field ir Error 


ERROR ATTEMPTING TO COMPUTE NEXT 
STORAGE LOCATION FOR NUMERIC LOCATION 
FIELD, BASIC SYMBOLIC ORIGIN HAS NOT 
BEEN SET. 


~— Field in Error 


UNDEFINED SYMBOL USED IN LOCATION 
FIELD. 


+— Field in Error 


MIS".NG COMMA IN MATRIX DEFINITION. 


—— Field in Error 
ILLEGAL CONVERSION CODE, 
—— Field in Error 


INVALID CiiARACTER IN VALUE FIELD. 


—~— Fieid in Error 


~— Card Image (80 Characters) -——= 


Type cf Input Error 
If the location subfield is not blank or the 
integer 1, the card is not stored; it is 
printed as shown, and the error 
termination procedure takes over, 


The first value input had a blank location 
subfield, and no symbolic origin had been 
established. 


_—— 


A matrix element is to be stored, but no 
location definition of the matrix has been 
made by a previous M or L card, 


— 


A matrix element is to be stored, but no 
size definition of the matrix has been 
made by a previous M or L card. 


— 


The first value input had a numeric 
location subfield. No symbolic origin had 
been established before this poirt. 


A symbol not previously defined was used 
in the location subfield. 


—_ 


During an attempt to define a matrix (M 
or L conversion code), the third character 
in the value subfield was not a comma, 


— 


An undefined conversion code was entered. 


—_ 


There is an invalid character inthe value 
field, 


—e 


i) 


SAMPLE GAIN1 INPUT FORM 


A.3 


A sample GAIL1 input form is shown below. 


6600 INPUT DATA 


PROGRAMMER 


iS 
orn 
rs 


PaGE 


OaTe 


-— « VERIFIES 


KEVPUNCHED 


VALUE 


LOCATION 


AUT TECAITT 


LOCATION 


atta 


AEROSPACE FORM 3770 REV §-67 


A-14 
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APPENDIX B 


TRACE DATA DECK STRUCTURE EXAMPLES 


B.1 ORBIT DETERMINATION RUN (ITIN=2) 
B.1i.1 Single- Vehicle with OBSERVATION Card Input 


The deck setup for an orbit determination run for one vehicle is shown in 


Fig. B-1. All observational measurements are input in the OBSERVATION 
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Fig. B-1. Teck Setup for a Single-Vehicle Orbit Determination Run 
with OBSERVATION Card Input 
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BiA.2 Multiple-Arc with OBSERVATION Card Input 


The deck setup for an orbit determination of more than one vehicle in which 
measurements from the many vehicles are used (but each observation is 
associated with only one vehicle) is shown in Fig. B-c. In this case, 

only one model group (MODEL, MASS, STATION, SENSOR, ATA, and 
CONSTRAINT data blocks) is input; it 1s followed by the VEHICLE blocks 
for ali vehicles, and then by pairs of REJECT and OBSERVATION blocks 


for each vehicle. 


Each OBSERVATION card should have a vehicle number (Sec. 15), and 
VEHID and BTIME are required in each VEHICLE block. For convenience, { 
in the VEHICLE block, the input data not overridden carries over from one 


vehicle to the next (characteristic of GAILi1 input). 


| 
I 


END QF INPUT 
TY 
OBSERVATION / 


¥ 
gs/ 
ee OBSERVATION 


7 


ae 
STATION 4 


I 5 


i 
Fig. B-2. Deck Setup for a Multiple- Arc Orbit Determination Run 
with OBSERVATION Card Input 
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Bei3 Single-Vehicle with Card Image Observation File Input \ 
The deck setup for an orbit determination run for one vehicle, using 
observational measurements found on a card image observation file (TAPE4) 
is shown in Fig. B-3. This setup is indicated by the input item BCDIN#0 
(Sec. 11.2.1) in the VEHICLE data block and the fact that the OBSERVATION 
card is followed by the END QF INPUT card. 
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Fig. B-3. Deck Setup fcr a Single- Vehicle Orbit Determination Run 
with Card Image Observation File Input 


B,1.4 Multiple-Arc with Card Image Observation File Input 


In an orbit determination run with observational measurements from several } 
vehicles (each item of data is associated with only one vehicle), the data deck i 
requires only one set of model group data blocks (Fig. B-4). For each 

vehicle there must be a VEHICLE block followed by the pairs of REJECT and 
OBSERVATION blocks. 


In the VEHICLE block, input data not overridden will carry over from one 
vehicle to the next. BTIME is required in each VEHICLE block to show the 
last integration time (more than likely, the last observation time). VEHID, 
which should match the vehicle number on the observation card images for 


that vehicle, is also required. 


For this case, observations are input by the card image observation file 
(TAPE4). This is indicated by BCDIN#0 (Sec. 11.2.1) and by a single 
zs @BSERVATION card for each vehicle. 
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3 | Fig. B-4. Deck Setup for a Multiple- Arc Orbit Determination Run 
% | with Card Image Observation File Input 
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B.1.5 Single- Vehicle Differential Correction with Binary 
Observation File Input 
Orbit determination for one vehicle with an input binary observation file 
(TAPE3) requires one meadel group of input data blocks and the VEHICLE 
block only from the vehicle group (Fig. B-5). The input item BTIME 
(Sec. 11.2.1) in the VEHICLE block contains the last observation time (or 
a later time) on the file. This time is considered in determining the final 
integration time. Note that no REJECT or OBSERVATION data blocks are 


allowed when a binary observation file is used. 
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Fig. B-5. Deck Setup for a Single- Vehicle Differential Correction 
Run with Binary Observation File Input 


e 
B.1.6 Multiple-Arc with Binary Observation File Input \. 
The deck setup for a multiple-arc orbit determination run with binary file 
observation input is shown in Fig. B-6. Observational measure.nents from 
many vehicles are used, but each item of data is associated with only one 
vehicle. One model group is input, and one VEHICLE block per vehicle is 
input. Each VEHICLE block contains tke input item BTIME (Sec. 11.2.1), 
which indicates the last observation time (or a later time) on the binary 
observation file (TAPE3) for this vehicle. These times are considered in 


determining the final integration time for each vehicle. 


Since all vehicles are assumed to be :ndependent, a complete set of appro- 
priate VEHICLE data should be provided for each. Input data not overridden 
is carried from one vehicle to the next. Again, note that no REJECT or 
OBSERVATION data blocks are allowed when a binary observation file is 


used. 
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5.2 EPHEMERIS GENERATION RUN (ITIN=3) W) 


B.2.1 Single- Vehicle 

The deck setup for a single-vehicle ephemeris generation run (Fig. B-7) 
consists of the following: the MODEL and the (optional) MASS data blocks 
from the model group and the VEHICLE block. The VEHICLE block must 
contain the input item PTIM (Sec. 11.3.1.1), which indicates the print times 


and the final integration time. 
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Fig. B-7. Deck Setup for a Single- Vehicle Ephemeris Generation Run 
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B.2.2 Multiple-Arc 


The deck setup for an ephemeris generation run with several vehicles 

(Fig. B-8) consists of one set of MODEL and the (optional) MASS data 

blocks and one block of VEHICLE data per vehicle. The input item PTIM 
(Sec. 11. 3.1.1) is specified in each VEHICLE block to indicate the print times 
and the final integration time for each vehicle. In the VEHICLE block, input 


data not overridden is carried over from one vehicle to the next. 
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Fig. B-8. Deck Setup for a Multiple- Arc Ephemeris Generation Run 


B,3 MEASUREMENT DATA GENERATION RUN (ITIN=4) a, 
B.3.1 Single- Vehicle 


The deck setup for a singie-vehicle measurement data generation run 
(ITIN=4) requires one set each of the model and vehicle groups of input data 
blocks, including MODEL, STATION, VEHICLE, DATA GENERATION, and 
the optional MASS and SENSOR blocks (Fig. B-9). 


The DATA GENERATION input block includes both DATA GENERATION I 
and II cards (Sec. 12) unles3 the VEHICLE block contains the input item 
JRIST=1 (Sec. 11.4.1). In this case, the DATA GENERATION II cards 


are omitted. 
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-ig. B-9. Deck Setup for a Single- Vehicle Data Generation Run 


B.3.2 Multiple-Arc 

The multiple-arc data generation deck (Fig. B-10) contains one set of model 
group data blocks and one set per vehicle of the vehicle group data blocks 
(VEHICLE and DATA GENERATION cards). If the VEHICLE block contains 
the input item JRIST=1 (Sec. i1.4.1), the DATA GENERATION II cards 

are omitted from the DATA GENERATION data block (Sec. 12). 


Note that all vehicles are assumed to be independent; therefore, a complete 
set of VEHICLE and DATA GENERATION blocks should be provided for each 
vehicle. For convenience, in the VEHICLE block, input data not overridden 


is carried over from one vehicle to the next. 
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Fig. B-10. Deck Setup for a Multiple- Arc Data Generation Run 
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B.4 COVARIANCE ANALYSIS RUN (ITIN=§} 


Deck setups for single-vehicle covariance analysis «uns always contain the 
same basic model groups of input data blocks (MODEL, MASS, STATION, 
SENSOR, ATA, and COVQ), whereas the vehicle data blocks (VEHICLE 

and DATA GENERATION) and observation data blocks (REJECT and 
OBSERVATION) are variable. Various deck setups for covariance runs are 


shown in Figs. B-11 through B-i5. 


Input/output options are specified in the MODEL data block by input items 
PPBPX and PRC@V (Sec. 2.5.1). The model P/Q parameter specifications 
are found in the input items PRAM, MPRAM, and GPRAM (Sec. 2.1.5). 


A priori input for the normal matrix Ata or the variance-covariance matrix 
(ras is specified in the ATA data block (Sec. 6). A priori input for 
Q-parameters is specified via tae COVQ data block (Sec. 8). 


Vehicle-depende..i P- and Q-parameters and the print schedule are specified 
by input items VPRAM (Sec. 11.1.14) and PTIM (Sec. 11.5.1). 


Observational measurement times may be specified by: 


OBSERVATION cards 

Card image observation file (TAPE4) 
Binary observation file (TAPE3) 
DATA GENERATION data block input. 


When the observational measurement times are specified by card image tape 
or card images, measurement times may be rejected by inputs to the REJECT 
data block. 


: T Toga t : 
In case of an update to an input A Aor (A A) = matrix, no measurement 


times are required [®PB@X(D) and PPBPX(E), Sec. 2.5.1]. 
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Fig. B-11. Deck Setup for a Covariance Analysis Run with 
OBSERVATION Card Input 
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Fig. B.12. Deck Setup for a Covariance Analysis Run with Card Image 
Observation File Input 
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Fig. B-13. Deck Setup for a Covariance Analysis Run with Binary 
Observation File Input 
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Fig. B-14. Deck Setup for a Covariance Analysis Run with Data Simulation 
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Fig. B-15. Deck Setup for a Covariance Analysis Update Cnly Run 


B.5 MULTIPLE ITINERARY 


B.5.1 Ephemeris and Data Generation (ITIN=34) 


The deck setup for a case in which an ephemeris generation and a data 
generation are desired for the same vehicle, using the same model group 
blocks, is shown in Fig. B-16. Included in the MODEL and VEHICLE data 
blocks are all input items required for both the ephemeris generation and 
the data generation runs, Even though an ephemeris generation run does 
not require STATION, SENSOR, or DATA GENERATION input blocks, they 
are included in the deck setup because a data generation run does require 


them. 
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Fig. B-16. Deck Setup for a Multiple-Itinerary Run (ITIN = 34) 


B.5.2 Multiple Itinerary (ITIN=3452345) 


The deck setup for Example B (Sec. 1.4.1) is shown in F.g. B-17. In this 
case, the initial ITIN functions 3, 4, and 5 cause the generation of the nomi- 
nal ephemeris, the look angles, and the covariance analysis, respectively. 
Reconstruction takes flace starting from the nominal initial values. When 
the iterations terminate, the trajectory for the converged solution is used to 


repeat the three processing functions. 


The MODEL and VEHICLE data blocks contain the necessary input items for 
all functions. The STATION cards are necessary for ITIN functions 2, 4, 
and 5, The optional ATA data block can contain Ata or (atay"! (Sec. 6). 
It is used for the orbit determination and/or the covariance analysis func- 
tions. The DATA GENERATION cards are necessary to generate the look 
angles, and the OBSERVATION block contains the measurements used in the 


covariance analvsis and orbit determination steps. 
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Fig. B-17. Deck Setup for a Multiple-Itinerary Run (I[TIN = 3452345) 
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B.5.3 Multiple Itinerary (ITIN=323) 
The deck sctup for Example A (Sec. 1.4.1) is shown in Fig. B-18. A 


trajectory is generated from the input initial conditions, and an ephemeris 
is generated and output. Then, another trajectory is reconstructed from the 
observational data, and comparable ephemeris output is printed. The 
VEHICLE data specifies the nominal initial conditions for the trajectory, 

the amount of printed output, and the vehicle parameters to be differ- 


entially corrected in the reconstruction. 


The observational measurements can be input by cards, card image file, or 
binary observation file (Sec, 11.2.1). Ifa card image file is used, the 
OBSERVATION block consists of a single card, with DBSERVATION in 
Columns 1 through 11. If the binary observation file is used, no REJECT 
or OBSERVATION data blocks are allowed. 
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Fig. B-18. Deck Setup for a Multiple-Itinerary Run (ITIN = 323) 
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B.6 STACKED CASES (ITIN=34 FOLLOWED BY ITIN=2) 
This data deck setup is used for Cases A and B of Appendix C (Fig. B-19). 


In Case A, an ephemeris is generated, and observational measurements 

are simulated. These measurements are written on a card image observa- 
tion file (TAPE4) and are then used in the orbit determination run of Case B. 
Note that in the first case, from the MODEL block to the first END OF 
INPUT, the data block setup is as described in Secs. B.2.1 and B.3.1. The 
VEHICLE data also includes input item ETAPE#0 (Sec. 11.4.1), so the 


observation2} measurements are wriiten on the card image observation file. 


The ITIN=2 function is consid*red an entirely new case because nothing is 

retained from the first case; new model, vehicle, and observation group 

inputs are therefcre included. The STATION input must be repeated, and 

the PBSERVA”’ION card must be followed by the END @F INPUT card. 

Input item BCDIN#0 (Sec. 11.2.1), specified in the second VEHICLE block, 

indicates that the observations are on the card image observation file, and 

BTIME (Sec. 11.2.1) indicates the time of the last observation. eee 
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Fig. B-19. Deck Setup for Stacked Cases (ITIN = 34, ITIN = 2) 
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B.7 SIMULTANEOUS VEHICLES 
B.7.1 Orbit Determination (ITIN=2, MULTV#0) 


The deck setup is shown (Fig. B-20) for many vehicles in simultaneous orbit 
that have the same model. The VEHICLE data blocks are input sequentially. 
Each requires a VEHID (Sec. 11.1.2) and the final integration time specified 
by BTIME (Sec. 11.2.1), but the integration times of the first vehicle must 
span those of any subsequent vehicles. Since the REJECT cards are time- 


dependent, not vehicle-dependent, only one block is input. 


Only one OBSERVATION block, containing the measurements from all vehi- 
cles, can be input. In this case, actual OBSERVATION cards are shown. 

If this data were input by the card image observation file (TAPE4), the block 
of cards would be replaced by a single OBSERVATION card. If it were input 
by the binary observation file (TAPE3), there would be no REJECT or 
OBSERVATION data blocks. 


The optional ATA block may contain an Ala or tata! matrix (Sec. 6). 
If MULTV=z2 or 3, the STAGE data block (Sec. 14) is optional, and if 
MULT V=2, the DEWM data block (Sec. 7) is optional. 
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Fig. B-20. Deck Setup for a Simultaneous-Vehicle Differential 
Correction Run with OBSERVATION Card Input 
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B.7.2 Data Generation (ITIN=4, MULTV40) i) 


The deck setup for generating data for vehicles simultaneouslv in orbit that 
have the same model is illustrated in Fig. B-21. Each VEHICLE data block 
contains a VEHID (Sec. 11.1.2). 


Only one DATA GENERATION block is input; it contains DATA GENERATION I 
(Sec. 12.1) ane II (Sec. 12.2.2) cards. 
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Fig. B-21. Deck Setup for a Simultaneous- Vehicle Data Generation Run 


5.7.3 Covariance Ant!ysis (ITIN=5, MULTV=1, 2) 


The deck setup for a simultaneous-vehicle covariance analysis run 
(Fig. B-22) is basically the same as that for a single vehicle (Sec. B.4). 
The main difference is that, since the vehicles are in simultaneous orbit, 


all VEHICLE data blocks are input before the DATA GENERATION, REJECT, 
or OBSERVATION data blocks. 


Each vehicle requires a VEHID (Sec. 11.1.2), and the first requires the 
input of PTIM (Sec. 11.5,1). The integration times of the first vehicle 


must span those of any subsequent vehicles. 
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C. SAMPLE OUTPUT DESCRIPTIONS 


INTRODUCTION 


TEST CASE A: ECI SINGLE-VEHIC LE EPHEMERIS GENERATION 
AND SIMULATED MEASUREMENT DATA GENERATION RUN 
(ITIN = 34) ........ 


TEST CASE B: ECI SINGLE- VEHICLE ORBIT DETERMINATION 
f RUN (ITIN = 2) 


TEST CASE C: ECI SIMULTANEOUS- VEHICLE COVARIANCE 
ANALYSIS RUN (ITIN = 5) 


TEST CASE D: MCI SINGLE- VEHICLE EPHEMERIS GENERATION 
RUN (ITIN = 3) ..... 


TEST CASE E: ECI SINGLE- VEHICLE, POWERED FLIGHT 
EPHEMERIS GENERATION RUN (ITIN = 3) .. 2.22.22 eee eee eceee 
TEST CASE F: ECI SIMULTANEOUS- VEHICLE ORBIT 

DETERMINATION RUN (ITIN = 2) 


TEST CASE G: ECI SIMULTANEOUS- VEHICLE DATA 
GENERATION RUN (ITIN = 4) 
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APPENDIX C 


SAMPLE OUTPUT DESCRIPTIONS 


C.1 INTRODUCTION 


The printed output produced by the TRACE Program is described here for 
typical differential correction, ephemeris generation, data generation, and 
covariance analysis runs. The samples that follow are of actual output 
listings, annotated to reference specific portions of the output data. Wherever 
applicable, the output items described are cross-referenced to the corres- 
ponding input definition section in this document. If an item occurs in more 


than one sample listing, a citation is given only for the first appearance, 


The sample test cases included here are: 


e Test Case A ECI single-vehicle ephemeris generation 
and simulated measurement data genera- 
tion run (ITIN = 34) 


e Test Case B ECI single-vehicle orbit determination 
run (ITIN = 2) 


© Test Case C ECI simultaneous-vehicie covariance 
analysis run (ITIN = 5) 


e Test Case D MCI single-vehicle ephemeris generation 
run (ITIN = 3) 


e Test Case E ECI single-vehicle, powered flight 
ephemeris generation run (ITIN = 3). 


e Test Case F ECI simultaneous-vehicle orbit determina- 
tion run (ITIN = 2) 


e Test Case G ECI simultaneous-vehicle data generation 
run (ITIN = 4) 
The symbol —~—~—— ,_used throughout the sample output, indicates that some 
output has been omitted. Note that many output pages without descriptive inter- 
ruption have been included only to provide continuity. Cases A through D were 


run on TRACE Version 6.96 and Cases E through G on Version 7.27. 


34) 


: ECISINGLE-VEHICLE EPHEMERIS 
GENERATION AND SIMULATED MEASUREMENT 
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APPENDIX D 


TRACE INPUT FORMS 


INTRODUCTION 


The engineering specification forms and load sheets that follow are very 
useful in expediting the preparation of TRACE input. Some of these forms 
bear former TRACE designations (TRACE-66 and TRACE66), which will 


be removed when the forms are updated. 


Ln ot ne ee ee 


as 


TRAJECTORY ANALYSIS AND PROGRAMMING DEPARTMENT 


TRACE66 (A0104a) ENGINEERING REQUEST FORM 


GENERAL INFORMATION: 


REQUESTER TELEPHONE ___e_ss«OATE 


PROBLEM NUMBER OR J O. QUE OATE 


PROJECT | 


TYPE OF TRACE RUN: 


EPHEMERIS GENERATION 
ORBIT DETERMINATION 
DATA GENERATION 
COVARIANCE ANALYSIS 


~AaAS 
wer 


GENERAL INSTRUCTIONS: 


1F A RERUN, JSE OATA DECK LABELEO 


NEW® OATA OECK TO BE SAVED: €) YES (} NO 


DISPOSITION (TO 8& COMPLETED BY PRODUCTION COORDINATOR) 


SET UP BY __ RUN 10 DATE 
TAPE SAVEO IF REQUIRED 


JOB DISPOSITION ANO RESULTS: 


NOTE: ALL SPECIFIED INPUT UNITS MUST BE CONSISTENT. 


TIME OF RUN: 
3 cp PP 
OATE JOB COMPLETED 


OATA OEFCK LABELED AS 


AEROSPACE FORM 4138 ALY ET) 


£ 
& 


oR 


Pa TEES ES TATE OI 


at asi eRet o> ecirorny mars oe athe SEEM FONE IP Et IT te ee RE 
oe 
- | TRACE 66 
x TF 
ORBIT AND FORCE MODEL SPECIFICATIONS 
RUN TITLE: 
INITIAL CONDITIONS: : 
EPOCH: INDICATE THE TYPE USED VALUES 
YEAR —— 
MONTH a 
Oay 
TIME ZONE WO LL 
HOUR eae 
MINUTE iF 
SECOND 
NOTE: UR IC TYP IS NEGATIVE, FHEEN THE ENIES OF THT INITIAL 
CONDITIONS ARE ASSEMED TO BE INTERNAL. 
LUNAR INITIAL CONDETIONS ARE TYPES TE THROUGH 18, CORRESPONDING 
TO THE ABOVE FARTIERITE RENCED TYPES T THROUGH 8, 
GRAVITY MODEL: CHECK ONE 
SPHERICAL () STANDARD 6, 6 () STANDARD 7. 6 () 
Cy. Cy. () STANDARD 8, 8 () OTHER (ATTACH) = () 
SUN-MOON ATTRACTIONS: No () YES (3 
DRAG COEFFICIENT: Cy 2 Ls 
SEGMENTED DRAG COEFFICIENTS: {Cohwes - 
{C)A/®), “ ‘7 MME (gy = ee Bg = ee te! 
(Cya/W), = tgs ___ MME CU SAO gc gg es 
(CHAD, es _ a" MME (pA We = - — ‘> ———— MME 
; (CA), = (ps see MME (yA Wy = oon 
oa C5 TABLES: INDICATE THE OPTION TO BE USED AND ATTACH TASLES IF ANY 
(.) NONE Od cae me () cy = 
(LINEAR WEIGHT LOSS: wo () yes () \ At vain = LS 
vw: LBS/MIN we _MMF iF MME 
INSTANTANEOUS WEIGHT CHANGES: NO ( ) YES \ ) AITACH TABLES 
ATMOSPHERE MODEL: = CHE( K ONE AND PPOVIDE INPUTS AS REQUIRED 
AKDC 1989 () US STD 1962 () 
I XONENTEA.. (AN Bie Ps _ 8 ee 5 — 
N‘I Lb Pan 1 al Se Oey = cy = ae 
oor eas an ie eee 
oe , ee 
LOCKE ED: JACCHIA ny O2 = FLUX = 
EMSe 1967 On Fy ae ee ee =o Swe ee 
JACCHIA 1904 (UNG rraies oe ce ene es (EDs lin ne . . 
JAC CHIEA 1964 (MOD) @), Fiat Fier” op 2 — 
RI. Ch Vio. Yor eee NG 
Wag yg g bp ORR ARF TOE EUSC TIONS OF HIME, AB EACH TAILS 
SOLAR RADIATION PRESSURE: no () yrs ©) 
s %, i, ® 
AW = HTL cy LBS/ET tae 
ONBOARD ACCELEROMETER: No () VES (0, ATTACH TABLES 
HINITE THRUSTING: No (} ves (0, ATFACH TABERS 
ORBIT ADJUSTS: no () VES © ATTACHE TARTES 
INPtY DENCIPGNITS ARE IN UU tem cm am cr etd 


AtRSP ACE POM 6134 REV 8-7) 


EASED MY, 


ea tahie 


aw 


TRACE66 
GRAVITY MODRL SPECIFICATION 


Check one: No normalization {) 
APL normalization ( ) 
KAUIA normalization ( ) 


— —— — —— —_. 
— —_— — ee —_ —- 
— —_ — —— —_ 
—— — - = — —w9a 
—_ ——— me _— = 
— —— — — AR 
— — — —— 
—— — ———————— oo =— —= ae: 
—o — —— — — — eer 
—_— — =e — —_ 
—o —_— — —_ — 
—o — — —= 
Yi 
—_— — — —— —— 
— — oe — — 
— — — — 
—— —_— — =— — — 
— — ee —_ — — 
—o —_— ——os — 
— — ——— — —— =— 
—_— — — —_—_ —_— ee a emt 
— — —— — — cee, 
— — — ——— 
— — —_ — 
— —_— —_— — 
-_—— _—_ — — 2 a eer 
—— —- —_ — 
— — — — 
— — —_—_— — — 


AEROSPACE FORM 4475 


cain iemminendianche dben es nh tainsnateaasdi Cine Reise et aa i a ia el Se a Ra RR ad a a 


TRACE66 
ATMOSPHERE AND Cy TABLES 


Check one: () ® = f(t); wath a maximum of $0 entries 
() kp = f(t), with a maxiwum of 50 entries 
Ck = f(c), with a maximum of $0 entries 


( ) Fuey = f(t), with a maximum of 12 entries 


() Cp = f(t), with a maximum of 25 entries 


(Cp = fh), with a maximum of 25 cotties 


(MME) of b(fe) f(t) ot f(h) UMME) of bf ft) f(t) of fh) 


AEROSPACE FORM 6669 


ONBOARD ACCELEROMETER: 


TRACE66 
EVENT SPECIFICATIONS 


T SCALE BIAS T SCALE BIAS 
MME FACTOR Fr/sec! MME FACTOR Fr/sec® 
FINiTE THRUSTING: 
TYPE P P, P P 
n y (FT/SEC) T, (I/MIN) NOT USED Vv (2 FT/SEC) 
2 Te C FT/SEC A/Go Vv (2 FT/SEC) 
5 A (FT/SEC) « B 
4 A (FT/SEC*) 6, (DEG) 6, (DEG) V (2 FT/SEC) 
3 A (FT/sEC4 6 (DEG) 6, (DEG) Vv (*FT SEC) 
6 T (LB) y 4 m n 
7 T (LB) ey (DEG) a (DEG) tv? FT/SEC) 
8 T (LB) ey (DEG) 4, (DEG) IV(SFT/SEC) 
Ts (MME) T, (MME) TYPE Pp PB y BY 
INSTANTANEOUS WiIGHT CHANGES: 
T(MMF) AWLBS) T(MMEF) AWOL BS) T(MME) AW(1. BS) 
ORBIT ADJUSTS: 
T(MMF) TYPE P, Py ty T° PE Pp P, r, 
pee ee 1 AR(E TSEC) AT(ET/SEC) AC(FT/SEC) 
———_— ee ee 2 K (FT/SEC) 6, (DEG) 4, (DEG) 
SS eee b) B (DEG) A (DEG) V(FT SEC) 


aCROSPACE FO 414t KEY ST 
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nos Cras Saale 


Aiea 


t 
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NASER ER E 


4 J 
See 


TRACE-66 
PARAMETER SPECIFICATIONS 


PARAMETER INITIAL 
IDENTIFICATION VALUE * BOUND 


Be ee TT ee Ta TAS eat le 
TAHT AALTTTEAA TI 
AIGA | | 


* THIS VALUE FIELD IS USED ONLY FOR OBSERVATION BIASES AND SCALE FACTORS. 


AEPOSACE POM 4148 


D-7 


STANDARD 
DEVIATION 


TRACE66 
ORBIT DETERMINATION SPECIFICATIONS 


ALL PARAMETERS MUST BE SPECIFIED ON THE TRACE66 PARAMETER SPECIFICATION FORM. 
OBSERVATIONS ARE TO BE INPUT BY; 


(>) CARDS 

() CARD IMAGE TAPE 

(.) BINARY OBSERVATION TAPE 

iF CARDS OR A CARD IMAGE TAPF ARE BEING USED AND THE OBSERVATION CATA TO BE USED IS A SUBSET 
OF THE DATA PROVIDED, THE BEGINNING AND END OF THE INTERVAL TO BE USED ARE GIVEN BY: 


BEGINNING YEAR MONTH DAY HOUR MIN SEC 


END YEAR MONTH DAY HOUR _ MIN SEC 


STATION INPUTS: ng AND ng ARE THE INDICES OF REFRACTION USED TO COMPUTE RFERACTION CORRECTIONS FOR 
ELEVATION AND RANGE DATA A ‘4/"* INDICATES THE NONINAL VALUES ARE TO BE USED, AO" 
INDICATES NO CORRECTIGNS AND IF A NON-STANDARD INDEX OF REFRACTION 1S TO BE USED, IT 
NUST BE PROVIDED. 


STA LAT LONG HEIGHT INTERFEROMETER Ne M, 
NAME (DEG) (DEG) (FEET) P Q (312 2 107) (330 « 10%) 


ae 


HN TH 


OBSERVATION WEIGHTS (0's): 


DATA DATA DATA DATA 
STA TYPE o 


SUPPLY THE REQUIRED INPUTS FOR THE FOLLOWING If OTHER THAN NOMINAL. 


MAXIMUM NUMBER OF ITERATIONS 
RELATIVE CONVERGENCE TEST (.001) 
ABSOLUTE CONVERGENCE TFST (0) 
a FOR a-SIGMA EDITOR (100) 


OPTIONS: CHECK THE OPTIONS DESIRED AND SUPPLY THE REQUIRED INPUT 


(.) SPEED OF LIGHT CORRECTION 

(.) POSITIVE CORRECTION 

(.) NEGATIVE CORRECTION 

PRINT OBSERVATION RESIDUALS 

PRINT OBSERVATION PARTIAL DERIVATIVES 

PRINT ORDERED (CORRELATION MATRIX 

PRINT CORRELATION MATRIX AFTER LAST ITF RATION 

SUPPRESS CORRELATION MATRICES 

PRINT ATA AND (ATA)"! AFTER EACH ITERATION, OR LAST ITERATION ONLY 
PUNCH ATA AFTER LAST ITERATION 

PUNCH SOLUTION VECTOR 

PRINTER PLOT OF OPSERVATION RESIDUALS AFTER LAST ITERATION 
DATARATE IS SEC 

(.) GENERATE AND SAVE OBSERVATION RESIDUAL TAPE 


~wew ewe we wr ww 


aGRosPaced Pome 6137 ety $7) 
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TRACE -66 
TRAJECTORY OUTPUT SPECIFICATIONS 


OUTPUT TINES: TIMES ARE GIVEN AS NINUTES FROM EPOCH ( ) OR MIDNIGHT ( ) 


FROM TO EvERY MIN 


SPECIAL OUTPUT TIMES: (y) 


SPECIAL 
B=% ©) LaT. LONG, ALT. 
MIN, MAX. ALTITUOE ( ) 
SPECIAL LATITUOES ( ) 
SPECIAL LONGITUDES ( ) 
SPECIAL ALTITUDES ( ) 
OBSERVATION TIMES ( ) ae 
OUTPUT REQUESTED: (y) 
STANOARO ECLIPSING 
4) EVERY TIME 0) 2) STANOARD () 
b) NOOES ONLY “y b) OPTIONAL 0) 
ELEMENTS VARIATIONAL EQUATIONS 
8) EVERY TINE =) *) STANDARO 0) 
b) NOOES ONLY ©) b) OPTION B () 
©) OPTION C () 
GEOMAGNETIC LAT. ANDLONG. ( ) GROUNO TRACK TAPE ) 
SUN-MOON ANGLES DIFFERENCE TAPE 0) 
) STANOARO 
b) OPTION Y 
¢) OPTION 2 a) 


AYRE Se DT PE 


EER 


AEROSPACE FORM 6140 REV 6-66 


TRACE66 
COVARIANCE ANALYSIS SPECIFICATIONS 


ALL PARAMETERS MUST BE SPECIFIED ON THE TRACE66 PARAMETER SPECIFICATION FORM. 


OBSERVATION GENERATION SPECIFICATIONS MUST BE SUPPLIED ON THE TRACEG6 DATA 
GENERATION SPECIFICATIONS - SHEET 1. 


IF OBSERVATIONS AR TO BE INPUT RATHER THAN GENERATED, INDICATE THE SOURCE 


€) CARDS { ) CARD IMAGE TAPE {) BINARY OBSERVATION TAPE 
IF CARDS OR A CARD IMAGE TAPE ARE BEING USED AND THE OBSERVATION DATA TO BE USED IS A SUBSET OF 
THE DATA PROVIDED, THE BEGINNING AND END OF THE INTERVAL TO BE USED ARE GIVEN BY: 


BEGINNING YEAQ MONTH . DAY HOUR MIN SEC. 


END YEAR MONTII OAY NoUR MIN. SEC. 


STATION INPUTS: Ne AND fp ARE THE INDIC FS OF REFRAC TION USFO TO COMPUTE REFRACTION CORREC TIONS FOR 
ELEVATION AND RANGE DATA. A‘Y" INDICATES THE NOMINAD VALUES ARE TO 8h USED, A “O" 


INDICATES NO CORRECTIONS, ANO IF A NON-STANDARD INDEX IS TO KF USED, IT MUST BE PROVIDEO. 


STA LAT LONG MEIGHT INTERFEROMETER on 1p 
NAME (DEG) {DEG) (FEET) P Q (312 10°) (350 2 10") 


OBSERVATION WEIGHTS (0's): 
OATA 
STA TYPE Q 


HT I 


OUTPUT TIMES: TIMFS ARE GIVEN IN MINUTES FROM EPOCH ( ) OR MIDNIGHT ( ) 
FROM To TVERY MIN 


OUTPUT REQUESTED: (') P Q 


€ (Pg) = (ATAY? TE ACS USESED THD ENTIRE COVARIANCE 


TQUINOC TIAL 
MIL ASUREMENT COVARTIANC I 


PRINT AP AQ 
SAVE PLOT TAPE 


©y a 
© ARTESIAN () C3) MATRIN BIEL BE PRINTED TE A “DY ds USED, 
ORHIT PLANE CX <a) ONLY THE SQUARE ROOT OF TIF DIAGONAL 
SPIE RICAL () () WIL BE PRINTED 
VURMENY (ih pales fe) 
PY RIOD, A) PER {) {) 

Ca O64 

Cis eee) 


an 


INPUT A PRIORI UNCERTAINTIES P PARAMETER Q PARAMETER 


NO A PRIORI INFORMATION 
ATA MATRIX GIVEN 
(ATAS?! MATRIX GIVEN 
o's GIVEN 

C(Q) MATRIX GIVEN 


~ aan 

~~~ 
~ 
- 


a 
~~ 


AEROSPACE FORM 41 Sy REV 5-71 
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Separs 


TRACE66 
DATA GENERATION SPECIFICATION SHEET 1 


STATIONS 
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